THE BELL SYSTEM 

TECHNICAL JOURNAL 

volume xxxii JANUARY number 1 

Copyright, 19BS, American Telephone and Telegraph Company 

Surface Properties of Germanium 

By WALTER H. BRATTAIN* and JOHN BARDEENf 

(Manuscript received September 3, 1952) 

The contact potential (c.p.) and the change of contact potential with il- 
lumination (Ac.p.) L of several germanium surfaces have been measured. 
The reference electrode used was platinum. It was found that the c.p. could 
be cycled between two extremes about 0.5 volts apart by changing the gaseous 
ambient. Ozone or peroxide vapors gave the c.p. extreme corresponding to 
the largest dipole at the Ge surface. Vapors with OH radicals produced the 
other extreme. There is a one to one correlation between c.p. and (Ac.p.) L . 
For 12-ohm cm n-type Ge (Ac.p.) L was large and positive when the surface 
dipole was largest, decreased to zero and became slightly negative as the 
surface dipole decreased to its smallest value. The variation for 12-ohm 
cm p-lype Ge was just opposite as regards both sign and dependence on 
surface dipole. The surface recombination velocity was found to be inde- 
pendent of c.p. For a chemically prepared surface it was 50-70 cm/sec 
and 180-200 cm/ sec for n and p-type surfaces respectively. A theory is 
given that explains the results in terms of surface traps, N a per cm 2 donor- 
type traps near the conduction band and N per cm 2 acceptor-type traps 
near the fdled band. A quantitative fit with experiment is obtained with 
only one free parameter. The residts are direct evidence for the existence 
of a space charge layer at the free surface of a semiconductor. 

INTRODUCTION 

Every one is familiar with the fact that it is necessary to expend 
energy to remove an electron from a conducting solid. This energy is 
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called the work function. The work function is caused in part by a 
charged double layer or dipole at the solid surface. In metals this dipole 
extends over a distance of the order of 10~ 8 cm. In semiconductors how- 
ever part of the dipole extends into the semiconductor to a distance of 
the order of 10~ 8 to 10 -4 cm depending on the properties of the semi- 
conductor. This part of the surface dipole is called the space-charge 
layer. The rest of the surface dipole has approximately the same extent 
as in metals. 

In the space between any two conducting solids there is a contact 
potential caused by the difference between the work functions of the 
two surfaces. One can measure this contact potential by several meth- 
ods. We have used an adaptation of the well known method of Kelvin. 
If one has a reference electrode whose work function remains constant, 
then by measuring the c.p. between this electrode and another surface 
one can measure any changes in work function or total dipole of this 
second surface. 

The above method has been used to study the properties of the 
germanium surface in a gaseous ambient at atmospheric pressure. It 
has been found that the total dipole at the germanium surface can be 
changed by changing the ambient and further that by proper control 
of the ambient the surface can be cycled back and forth between two 
extremes of small or large dipole corresponding to a c.p. change or work 
function difference of the order of one-half volt. 

If one upsets the thermal equilibrium in the germanium by creating 
excess electron-hole pairs near the surface, the potential of the surface 
will change until a steady state is reached. When the extra electron- 
hole pairs are introduced by illuminating the surface with light, the 
potential change shows up as a measurable change in contact potential 
between the reference electrode and the Ge surface. 1 It has been found 
that this contact potential change on illumination (Ac.p.) t is large and 
positive on n-type Ge when the surface dipole is large, then decreases 
to zero and becomes slightly negative as the surface dipole decreases 
to the smaller extreme. For p-type Ge the (Ac.p.)ji is large and negative 
when the surface dipole is small and goes through zero and becomes 
slightly positive as the surface dipole increases to the larger extreme. 

One can describe qualitatively what is going on as follows. The extra 
hole and electron pairs created by the light diffuse either to the interior 
or to the surface to recombine. The recombination in the interior is 
governed by the body life time t. The surface recombination is charac- 
terized by a recombination velocity v, . When the surface is illuminated 
its potential, with respect to the interior, changes until the combined 
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flow of holes and electrons to the surface and interior is just equal to 
the rate they are being created by the light. The sign and magnitude 
of the potential change for a given illumination depends on the body 
properties of the germanium and on the size of the space charge layer. 

These experimental results are direct evidence for the existence of a 
space charge layer at the free surface of a semiconductor. They not only 
confirm the results obtained for silicon surfaces 1 but go much further 
in that they enable one to determine how the layer is changed by the 
gaseous ambient used. 

It is known that the surface recombination velocity, v„ , can be 
changed, by large factors, by surface treatment. 2 For mechanically 
treated surfaces v, approaches thermal velocities. Every hole or electron 
striking the surface recombines. For such a surface it is found that 
(Ac.p.)j, is too small to be measured. On the other hand v s can be as low 
as 100 cm/sec for chemically polished or etched surfaces such as those 
used in the experiments where (Ac.p.) L was measured. In this case one 
wishes to know how v a depends on the gaseous ambient. This was meas- 
ured for the same surface used in measuring (Ac.p.)r, and it was found 
that, for the ambients used, v, is approximately a constant and there- 
fore independent of the other surface changes. 

A quantitative theory, some details of which are in the Appendix, has 
been formulated to explain the results. It is proposed that there are 
two types of recombination traps at the surface : donor type, N a per cm 2 , 
with energies, E a , near the conduction band and acceptor type, Nb per 
cm 2 , with energies, Eb , near the filled band. Surface recombination 
takes place by electrons and holes successively going into one of the 
two types of traps. To account for the fact that v, is unchanged by 
changes in ambient, it is assumed that the concentrations of these traps 
are independent of ambient. Changes in c.p. with ambient are assumed 
to result from adsorption and desorption of fixed ions which are at an 
effective distance f~2X 10~ 7 cm outward from the surface traps. A sche- 
matic energy level diagram is given in Fig. 13, to be discussed later. 

The charge of the ions is compensated mainly by charges in the 
surface traps which, together with the ions, form a double layer. A 
large part of the change in c.p. with ambient results from changes in 
this double layer. There is also a change in barrier height, —eV B , 
associated with the redistribution of electrons in the traps. An increase 
in negative ions on the surface requires a decrease in number of electrons 
in traps, and thus a higher barrier. 

Part of the change with light, (Ac.p.) L , occurs in the body of the 
semiconductor and part occurs across the barrier layer. Changes in V B 
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and occupancy of the traps compensate in their effect on surface re- 
combination, so that v, is unchanged by ambient. This means that 
changes in concentrations of electrons and holes in the interior with 
illumination and thus the body contribution to (Ac.p.)z, are independent 
of ambient. Changes in (Ac.p.) L with ambient result soley from changes 
in V B and are in the directions we have described earlier. 

The concentrations of carriers, n and p, and their change with light 
are obtained as follows. Drift mobilities, /x» and \i p , and the equilibrium 
product, np, are known from earlier experiments of G. L. Pearson, J. R. 
Haynes and W. Shockley. 3 From these, and the resistivity of the sample, 
which was measured, n and p can be determined. The light source is 
calibrated in terms of hole-electron pairs created per cm per sec. 
The recombination rate is determined from the body life time r and the 
surface recombination velocity. From these latter three measurements, 
one can calculate the steady state density of electrons and holes near 
the surface when it is illuminated. 

Mention should be made here of the fact that oxygen has been found 
to play a definite role on the Ge surface. The large extreme in dipole is 
obtained when active oxygen (ozone) is introduced into the gaseous 
ambient. Peroxide vapors have the same effect. The other extreme is 
produced by vapors having an OH radical, water vapor, alcohol etc. A 
number of vapors not falling in either of the above classes have little 
or no effect on the surface dipole. Another result is that the difference 
in work function or dipole between n- and p-type Ge is small. This is 
to be expected from previous work. 

We shall first discuss the experimental technique and then give the 
experimental results. The main conclusions of the theory will then be 
outlined and compared with these results. 

EXPERIMENTAL METHOD 

The Ge surface to be measured and the reference electrode are mounted 
under a bell jar. Oxygen or nitrogen as desired is allowed to flow through 
the bell jar at a rate of approximately 2 liters per min. The volume of 
the bell jar is 16 liters. The gas used flows over a drying column of silica-gel 
and then calcium chloride. Means are provided for bubbling this gas through 
any desired liquid before it enters the bell jar. A spark discharge can be 
run in the gas flow line. The reference electrode is placed parallel to the Ge 
surface about 1 mm away. It is mounted on a vibrating reed which is 
driven, electromagnetically, at its resonant frequency of about 90 cycles 
per second and at an amplitude of the order of 0.1 mm. This varies the 
capacity sinusoidally giving rise to an electrical signal when any potential, 
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contact or other, exists between the surfaces. When the proper dc potential 
is applied between the surfaces this signal goes to zero. If no other poten- 
tials are present this dc potential is equal and opposite to the contact 
potential. A phase reference method* is employed to determine this balance 
point with a relative accuracy of ±5 X 10 -4 volts. A diagram of the Ge- 
reference electrode circuit is shown in Fig. 1 . Care must be taken to shield 
this circuit. Stray capacity reduces sensitivity and should be minimized. 
Charged insulators inside the shield will produce an apparent c.p. All con- 
ducting surfaces other than the Ge should be relatively far from the moving 
reference electrode. The surface is illuminated through a compound 
lens system by focusing the filament image, of a suitable projection 
bulb, on the germanium surface. This light passes through the grid of 
the reference electrode which removes about 10 per cent of the light. 
The light can be modulated by a square wave chopper, so that (Ac. p.)*, 
can be measured on an ac basis. Both Ta and Pt reference electrodes 
have been used. The Pt electrode appears to be somewhat more constant. 
If the Ge surface is replaced by a gold electrode the contact potential 
difference is practically independent of the changes in the gas ambient. 
The arrangement is such that two samples can be mounted in the bell 
jar and the reference electrode moved from one to the other without 
opening the bell jar. In this way two surfaces can be compared, without 
any question arising of long time drifts in the reference electrode. 

EXPERIMENTAL RESULTS 

/. Change of c.p. between Ge and Pt reference electrodes as a function 
of the gaseous ambient 

When this work was started the object was to find some means of 
varying the c.p. It was thought that the actual values of the c.p. would 
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Fig. 1 — Schematic of experimental circuit. 



H. R. Moore designed and made the electronic equipment used to do this. 
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be highly dependent on past history of the Ge sample. If however one 
could measure one or more other properties, such as (Ac.p.)z, and v 8 , 
on the same surface at the same time, then one could look for correla- 
tions between these properties. In this manner one might be able to 
eliminate the past history as a factor. From previous experiments in the 
highly variable ambient of room air we knew the order of magnitude 
of the variation to be expected. At first it was impossible to produce 
this range of variation under the bell jar. The contact potential always 
drifted in the direction of a positive extreme, i.e., small total dipole at 
the Ge surface. The only way found to get a really large change in the 
opposite direction was to lift the bell jar and expose the sample to room 
air. These phenomena were finally traced to the presence of negative ions, 
possible salt ions, in the room air. These were not present in the oxygen 
and nitrogen supplies we used for creating the ambient in the bell jar. 

It was found that the opposite c.p. extreme could be produced, under 
the bell jar, by running a spark discharge in dry oxygen as it was flowing 
into the system. The next step was to cycle the c.p. from one extreme 
to the other and back again. The procedure was to start with the spark 
discharge in dry oxygen, change to either wet 2 or wet N 2 and to end 
with dry 2 . The development of this dependable and reproducible 
cycle was a great aid to the proposed study. Fig. 2 is a plot of contact 
potential versus time for a single crystal slice D of Ge cut from a melt 
that was p-type. The surface was prepared by removing some of the 
Ge with a silicon carbide (180 mesh) blast of approximately ten pounds 
air pressure. The Ge was then mounted in the bell jar within one-half 
minute of the "sandblast," and the dry 2 flow started. The c.p. was 
followed for a few minutes to be sure everything was working properly. 
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Fig. 2 — Contact potential cycles for sandblasted sample D. 
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Fig. 3 — Contact potential cycles for etched sample A. 



In Fig. 2 zero time is taken after the spark discharge was run in the 2 
flow line for 2 minutes. This started the first cycle. After approximately 
17 minutes the O2 was made to bubble through H 2 0. Fifteen minutes or 
so later the 2 was changed back to dry. At this time, 32 minutes, the 
flow rate was increased by a factor of three. The c.p. was followed for 
about 17 minutes, then the process was repeated. The results and the 
reason for the choice of time intervals, are all evident from a study of 
Fig. 2. The spark discharge in 2 decreases the c.p. After this treatment 
the c.p. increases with time, most of the change occurring in the first 
15 minutes. The wet 2 then increases the c.p. to a maximum value 
which is reached in about 15 minutes. Finally the dry 2 reduces the 
c.p. It is evident that there is a quasi-equilibrium value of c.p. in dry 
O2, to which the c.p. returns after either extreme treatment. At first 
there is quite a large shift from cycle to cycle but this shift gradually 
disappears as the cycling is continued. In Fig. 2 cycles 2, 4, 6 and 10 
are shown. Very little change takes place after cycle 10. Such results 
have been obtained many times over a period of two years. When al- 
lowance is made for shifts in work function of the reference electrode 
and for the fact that the experimental technique improved as the work 
progressed it is found that all the results for a given sample of Ge are 
very consistent. 

In Fig. 3 are shown similar results for an n-type slice A. In this case 
the surface was first ground or sandblasted to remove any films, and 






8 THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1953 

0.8 

0.7 

f2o.e 



Z 0.5 



• 



£ 0.4 




0.3 



20 24 28 32 

TIME IN MINUTES 



Fig. 4 — Contact potential cycles for etched sample D. 

then given a polishing etch (CP-4). After the etch the surface was 
washed in running distilled water of reasonable quality. The surface was 
then dried with filter paper and kept covered until placed in the bell 
jar. 4 From here on the procedure was the same as before. In this case 
cycles 1, 3 and 11 are shown. The surprising result is the similarity be- 
tween Fig. 2 and Fig. 3. To a first approximation one is tempted to say 
that the dipole of a Ge surface, in the bell jar atmosphere is independent 
of past history. Within certain limits this is approximately true. There 
are differences between Fig. 2 and Fig. 3 but they are small and probably 
due to differences in surface treatment. Fig. 4 shows the results for 
p-type slice, D, when this surface is etched as above. This slice and the 
slice A were placed in the bell jar at the same time. Cycle 1, Fig. 3, was 
taken on A, cycle 2, Fig. 4, on D and so on. Any differences between 
the results in these two figures are to be attributed to the differences in 
samples. They cannot very well be ascribed to the reference electrode 
and the initial surface treatments were as nearly the same as they 
could be made with reasonable care. By making runs of this type two 
samples at a time, different samples and different surface treatments 
can be intercompared. This method eliminates the shifts in the work 
function of the reference electrode that sometimes occur from run to 
run. Such results can be illustrated by plotting the data for different 
samples and different surface treatments for cycles 10 or greater where 
the results for successive cycles are the same. This has been done in Fig. 
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5 where one cycle each is plotted for both sample A and D, for each of the 
two sin-face treatments, sandblasting and etching. The curves are approxi- 
mately all of the same shape so that the differences between them can be 
described by giving the shift in contact potential necessary to superimpose 
the curves. This treatment works very well except for the first part of the 
cycle after the spark coil where the shift necessary is sometimes more and 
sometimes less than that needed to make all of two curves superimpose well. 
Note in Fig. 5 that the contact potential for sample A etched is always 
greater than for sample D etched. In the case of the sandblasted surface 
just the reverse is true. Also the contact potential for the etched surface of 
either sample is always larger than for the sandblasted. 

Using these methods comparable results were obtained on two other 
?i-type samples, C and E, of increasingly lower specific resistance. 
Taking advantage of the relation between carrier concentration and 
the position of the Fermi level we have plotted in Fig. 6, the contact 
potential for both the etched and the sandblasted surfaces versus the 
position of the Fermi level (E F — Ei) in electron volts. The contact 
potential values used were taken from the saturation values in wet O2 
but the shapes of the curves would be much the same for any other 
point in the cycle. The contact potentials are thought to be accurate 
to about 0.01 volts. Solid lines have been drawn through the points. 
Note that the contact potential for the etched surface is always greater 
than for the sandblasted surface, i.e., the work function is always less 
and that this difference is greater when E F = Ei or the germanium is 
nearly intrinsic. 
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Fig. 5 — Comparison of cycles for samples A and D with sandblasted and 
etched surfaces. 
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Fig. 6 — Contact potential dependence on position of the fermi level. 

For variety Fig. 7 shows results on c.p. for a sandblasted polycrystal- 
line sample of silicon. It is apparent that similar phenomena are taking 
place on the silicon surface, however the behavior is different. These 
preliminary results are shown to illustrate the generality of this method 
for investigating semiconductor surfaces in a controlled gaseous ambient. 

The time between cycles was somewhat variable. The first and second 
cycles were always taken immediately after the specimen was placed 
in the bell jar. After this, successive cycles were taken one or two a 
day over a period of a week or more. The bell jar was not opened during 
a run and a small flow of dry gas was maintained between cycles. Little 
if anything occurred during these idle periods, showing that the changes 
that did take place were due to the cycling. 

II. Change of contact potential with illumination 

When the Ge surface is sandblasted the change of contact potential 
with illumination (Ac.p.)t is too small to be measured. If however the 
surface has been prepared by the polishing etch, the (Ac.p.) L is easily 
observed. This change if not too small can be measured by finding the 
balance on the potentiometer for light off and light on. When the con- 
tact potential is changing with time, or if the change is small this is 
difficult to do. A better procedure is to chop the light at a definite 
frequency and measure the amplified output on an ac meter. This gives 
a continuous reading that can be read easily at any given time. If a 
filter is used to pass only those frequencies near the fundamental of 
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the chopping frequency, the improvement in signal to noise enables one 
to read very small changes. 

A plot of (Ac.p.)t in volts versus the contact potential for samples A 
n-type and D p-type is shown in Fig. 8. For the n-type sample the signal 
is large and positive when the contact potential is small and becomes 
small and negative as the contact potential increases. Except for the 
shift in the contact potential where (Ac.p.) L goes through zero, the re- 
sults for the p-type sample are practically the opposite of those for the 
n-type sample. Similar curves have been obtained cycle after cycle in 
many complete runs. Within experimental error the curves have the 
same shape for a given sample for all cycles in all runs. Sometimes the 
curve for the first cycle in a run will differ in shape from the rest. How- 
ever there are shifts such that the c.p. for which the light goes through 
zero (c.p.)o does vary from cycle to cycle throughout a run. Fig. 9 is a 
plot of (c.p.)o versus cycle number for p-type sample D and n-type 
sample A. The data are plotted for two distinct runs. Both of these 
units were in the bell jar when the measurements were taken. Consistent 
results of this kind give one confidence that the Pt reference electrode 
is staying constant. 

These experimental results can be summarized as follows. All data 
for a given sample can be superimposed by shifts in contact potential 
scales for the different cycles. A plot of (Ac.p.)r, versus c.p. — (c.p.) 
can be represented by a single curve for each unit. Moreover all the 
curves for all units both n- and p-type have quite similar shapes provided 
in the latter we plot [— (Ac.p.)xJ versus (c.p.)o — c.p. The shape of this 
curve is shown in Fig. 8 for a given sample and Fig. 9 shows how (c.p.)o 
varies from cycle to cycle. These plots adequately describe the results. 
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Fig. 7 — Contact potential cycles for silicon. 
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Fig. 8 — Change of contact potential with illumination versus contact po- 
tential, samples A and D etched. 

Not shown in Fig. 8 because of the scale used is the result that (Ac.p.)t 
for n-type material does not increase indefinitely as c.p. decreases but 
approaches a maximum value. Likewise — (Ac.p.)z, for p-type approaches 
a maximum as c.p. increases. Figures illustrating these results will be 
discussed after the theory is presented. 

Some of the experimental details require discussion. It is necessary 
to calibrate the ac response in terms of absolute potential change. This 
can be done by comparing the ac reading with the dc reading when the 
light signal is large. One can also do this by introducing a known square- 
wave signal across the potentiometer in Fig. 1, and reading the ac signal 
out. Both methods agree when allowance is made for variation of the 
light signal with frequency. The latter response is almost flat from 25 
to 300 cycles, but there is evidence for some very low frequency com- 
ponents in the dc measurements of (Ac.p.)r, . "When the light signal 
goes through zero the signal is small and the dc value is difficult to read. 
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At times evidence has been obtained to indicate that the dc value changed 
sign. At other times the dc (Ac. p.)*, behaved as if the place where it 
went through zero was shifted in c.p. from the point where the ac 
signal goes through zero. In view of this it was necessary to prove that 
the ac signal was changing phase at the zero point and not just going 
down in the noise and then increasing again without phase change. This 
was done by comparing the phase of the signal with a signal from a 
photocell placed in the same chopped light beam. By this means it was 
proved conclusively that the ac light signal was actually going through 
zero. 

Some data were obtained on change of contact potential with light 
on n-type samples C and E having progressively smaller specific re- 
sistances. It was found that (Ac.p.)z, decreased with specific resistance. 
It also decreased into the noise as the contact potential was increased, 
so that it could not be determined if it changed sign as for samples A 
and D. Because of the smaller signal (Ac.p.)z, could not be measured 
easily except by the ac method and so far the signal has not been cali- 
brated properly. 

Some preliminary data on p-type silicon indicate that (Ac.p.)*, for 
this sample was negative and that it decreased as c.p. was decreased. 
The magnitude of (Ac.p.) L for the same light intensity was much larger 
than for germanium and so far it has not been found to go through zero 
and change sign in the experimental range. 

III. Other methods of varying the c.p. 

In some cases N 2 was used in place of 2 . A spark discharge in the 
N 2 had very little effect on the c.p. On the other hand wet N 2 produced 
much the same effect as wet 2 . The positive extreme in c.p. was about 
0.1 volt greater in the case of wet N 2 . After the wet treatment dry N 2 
was not nearly as effective as dry 2 in reducing the c.p. to its inter- 
mediate value. This can hardly be due to a difference in dryness of the 
two gases since the same drying column was used in both cases. The 
results indicate that dry N 2 tended to leave the surface in whatever 
condition obtained before the dry N 2 flow was started, and that 2 
counteracts the effect of H 2 0. 

With dry N 2 as a carrier, other vapors were tried. A. N. Holden sug- 
gested trying a peroxide and picked out ditertiary butyl peroxide as 
being reasonably safe. Use of this vapor was found to produce the same 
changes as the spark coil in the 2 flow. Other vapors having OH radicals 
such as methyl alcohol and acetic acid were found to act the same way 
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Fig. 9 — Contact potential for zero light effect, (c.p.)o , versus cycle number, 
two runs for each sample. 

as water vapor. To prove that this was not caused by traces of water in 
the alcohol or acid the N 2 was bubbled through water solutions of 
H 2 S0 4 . These results indicated that one needs appreciable amounts of 
water vapor to produce the effect, much more than could be present in 
the alcohol or acid. Other vapors, such as carbon tetrachloride, methyl- 
chloride, nitrobenzene and ether, were found to have no effect on either 
the contact potential or (Ac.p.) L . Acetone has a small effect in the same 
direction as water. This is to be expected because this compound exists 
in part in a tautomeric form having an OH group. Vapor from 30 per 
cent H 2 2 , 70 per cent H 2 acted at first like a peroxide vapor and 
with a longer time of exposure behaved like water vapor. Small amounts 
of Cl 2 gas in N 2 produced the same change as the spark discharge in 2 
and after 14 minutes of flushing the bell jar with N 2 produced an ad- 
ditional effect when water vapor was introduced. On w-type samples 
before the usual increase in c.p. and decrease in (Ac.p.)z, there was a 
large increase in (Ac.p.) L . This was attributed to the reaction between 
the water vapor and the CI left on the Ge surface, producing oxygen. 
The nature of the change was a rapid shift in (Ac.p.)o and thus a mo- 
mentary increase in (Ac.p.)x; . This effect is only large in the first cycle. 
It indicates that the shifts in (Ac.p.)o plotted in Fig. 9 are probably due 
to an oxidation of the Ge surface as the cycling progresses. In all these 
experiments the relation between the c.p. and (Ac.p.)x, was essentially 
the same as that obtained in the standard cycle. 

One can detect the presence of thin surface films by electron diffrac- 
tion techniques. R. D. Heidenreich took electron diffraction pictures of 
a germanium surface immediately after the polishing etch and water 
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wash. He found the surface to be quite clean, with a film thickness less 
than 10 A. He also took pictures after the germanium surface had been 
cycled about fifteen times and found in this case definite evidence for a 
thin surface film. The film was either amorphorus or composed of very 
small crystals. He estimated the thickness to be between 20 and 50 A. 

In the discussion of the experimental work it has been assumed that 
all the changes in c.p. are to be attributed to the Ge surface and not to 
the Pt. It is not easy to give a definite proof that this is true. The fact 
that it is a reasonable assumption is suggested by the nature of the re- 
sults themselves. Almost identical results were obtained using a Ta 
electrode. In Fig. 10 we show the results of two cycles when the Ge was 
replaced with gold. While there are some changes they are almost an 
order of magnitude smaller than the changes when Ge is present. It 
follows that one would expect much the same results with either Ta, 
Pt or Au reference electrodes and that most of the changes are due to 
the Ge. No change of c.p. with illumination is observed when both 
electrodes are metals. In one case a small amount of H 2 was added to 
the N 2 flow. Here the c.p. between Pt and Ge changed rapidly but the 
(Ac.p.)/, did not change at all. Subsequent runs using the regular cycle 
indicated that the c.p. scale had been shifted corresponding to a reduc- 
tion in work function of the Pt. Except for this shift the results were 
the same and the shift disappeared in about one day. The conclusion 
was that H 2 had little effect on Ge but reacted with the Pt decreasing 
its work function. This is a good illustration of the power of this method 
of measuring more than one property of a semiconductor surface at the 
same time. If only c.p. had been measured the conclusions would not 
have been so clear cut. 

If one knew the work function of the Pt electrode then one would 
know the work function of the Ge. Work functions are all measured in 
high vacuum. We have not been able to think of a method of determin- 
ing the work function of any electrode in a gaseous ambient unam- 
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Fig. 10 — Contact potential cycles when germanium is replaced by gold. 
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Fig. 11 — Dependence of the velocity of recombination, i/,,on the contact 
potential cycles. 

biguously. It may well be impossible to do this in Bridgman's operational 
sense. Most physicists would agree that the Pt reference electrode 
probably has a work function of around 5 to 7 volts under these con- 
ditions, but this of course does not help much. 

IV. Measurement of surface recombination velocity v. 

At first we tried various methods of measuring v.'on the Ge surface 
at the same time that c.p. and (Ac.p.)i, were measured. No method 
was found that could be trusted. As the study progressed we realized 
that the c.p. of a Ge surface could be cycled in a reproducible way. 
Since the proper geometry for measuring v, was a rod or filament, 
rods were prepared of samples A and D with the same chemical surface 
treatment. The decay of hole electron pairs, created by a point light 
source, was measured as a function of distance from the light. If the 
body life time t and the dimensions of the filament are known one can 
then determine v„ . 2 These measurements were made while the gaseous 
ambient was cycled in the same maimer as before. The results for fila- 
ments cut from samples A and D are plotted against time in Fig. 11. 
The ambient atmosphere was changed as a function of time just as it 
was in Figs. 2 to 4. The main conclusion is that within the accuracy of 
this experiment there is no evident dependence of v, on the changes in 
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gas ambient and therefore no dependence on the corresponding changes 
in surface dipole. This result was somewhat unexpected and the first 
time the experiment was performed it was hard to believe that the 
previously measured changes in c.p. actually were taking place. In this 
case the gas ambient was experimented with to try to change v a and it 
was found that v, could be changed from the order of 10 2 cm/sec to 
greater than 10 5 cm/sec and back again by exposure of the filament to 
(NH) 4 OH fumes and then HC1 fumes respectively. 5 While interesting, 
this has no direct bearing on the other experiments. The experiment 
was performed again with freshly prepared rods and this time the 
cycling used in the c.p. experiments was rigidly adhered to, giving v s 
equal to 70 cm/sec and 200 cm/sec approximately for samples A and 
D respectively. The experiments were then repeated again using new 
filaments cut from the samples as close as possible to the surfaces used 
in the c.p. experiments leading to the results shown in Fig. 11, namely, 
v, equal to 50 cm/sec and 170 cm/sec respectively. From these experi- 
ments it was concluded that v, is approximately constant in the range 
involved and is determined by the nature of the sample and the surface 
treatment used. It was noted in some of these experiments that v, for 
the first cycle was somewhat larger than for the subsequent cycles. 
This change when it occurred is probably to be correlated with the 
changes in the early cycles in the c.p. measurements. Similar measure- 
ments on sample C gave v, equal to 1500 cm/sec. No measurements of 
v, were made on samples B and E. 

V. Other experimental measurements 

The specific resistance of each sample was measured near the surface 
used in the experiments. It was approximately constant across the 
surface but did vary slowly with depth in some of the samples. 

The body life times were measured on each sample. The thickness of 
the slices used was intentionally made large compared to their cor- 
responding diffusion lengths, about 0.5 cm for A, B, C and E and about 
2.0 cm for D. The mobilities were taken from J. R. Haynes 3 measure- 
ments: fi n = 3600 and n P = 1700 cm 2 / volt sec. There is some uncer- 
tainty as regards the exact value of the equilibrium product of holes 
and electrons, np, at 300°K. We have used the value 6.3 X 10 26 obtained 
from some unpublished data of G. L. Pearson. 

The light source was calibrated by replacing the germanium sample 
with one of F. S. Goucher's n-p junctions. 6 The bell jar, with every- 
thing else including the reference electrode, was left in their normal 
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positions. An average was taken over the filament image, and the ef- 
fective area of the p-n junction was determined and allowed for. This 
was done for all light intensities used. Most of the experiments were 
performed with a fixed intensity and the averaged result for this in- 
tensity was 6.0 X 10 15 , hole electron pairs per cm 2 sec. The rate of pair 
production was found to be proportional to the light intensity. 

Since practically all the light is absorbed in a depth (10~ 4 cm or less) 
that is small compared to the diffusion length, it is a simple matter to 
calculate the steady state increase 8p in hole electron pairs due to the 
light. The relation is 



8p = N/(v. + v d ) 



(1) 



where N is the rate of pair production, v a is the velocity of recombination 
at the surface and v d is the diffusion velocity for the minority carrier. 
Since N is proportional to light intensity it follows that Sp is too. 

The magnitude of (Ac.p.)t should depend on the light intensity. 
One might at first expect it to be proportional to light intensity. That 
this is not the case is shown by curve 1 in Fig. 12. Curve 1 is a plot of 
(Ac. p.)/, versus 8p for unit D on a log-log scale. A smooth curve has 
been drawn through the experimental points. As we shall see later, 
theory predicts that if (Ac.p.) L is large, it should be proportional to 
ln(l + 8p/a) where a is the equilibrium density of the minority carrier, 
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Fig. 12 — Dependence of contact potential change with illumination (Ac.p.) L 
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n in p-type material and p in w-type. That this prediction is borne out 
is shown by how well the experimental points fit the straight line curve 
2 in Fig. 12 where (Ac.p.)z, is plotted versus this quantity. The scale 
for 8p is shown along the bottom of Fig. 12 and that for tn{\ + 8p/a) 
along the top. Similar results were obtained for unit A. 

In Table I the parameters, specific resistance p in ohm cm, life time t 
in microseconds and the surface recombination velocity v, in cm/sec 
are given for each unit used. Also given are some pertinent quantities 
derived therefrom, namely the equilibrium densities of electrons and 
holes, n and p in number per cm and the increase in density at the 
surface 8p in number per cm 3 when the rate of pair production due to 
the illumination was 6.0 X 10 per cm sec. 

THEORY 

The constancy of v s throughout the range of surface dipole investi- 
gated puts rather stringent requirements on any theoretical model to 
be constructed. W. Shockley and W. T. Read have investigated the 
theory of recombination via traps. It is evident from their work that if 
one assumes a trap density peaked near a single energy and very small 
elsewhere, then v a will be constant over a range of surface dipole values 
provided that the peak energy is either near the conduction band or 
near the filled band. The experimental results make it appear very 
unlikely that the trapping mechanisms on the n and p-type surfaces 
are essentially different. It is assumed that both types of traps are 
present on the surface and that the traps are approximately the same 
for both n- and p-type samples. Further, it is assumed that the traps 
No of energy E a near the conduction band are donor type, i.e., neutral 
when filled and positively charged when empty. Likewise the traps N& 
of energy E& near the filled band are assumed to be accept or- type traps, 
i.e., neutral when empty and negatively charged when filled. The 
absolute charge on the traps is not important, however, because we are 
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concerned only with differences between the occupied and unoccupied 
states. 

The fraction of the traps which are occupied depends on the trap 
energy and on the position of the Fermi level at the surface. The latter 
in turn depends on the condition that the surface as a whole be neutral. 
A very obvious mechanism for changing the total surface dipole is the 
adsorption or desorption of ions on the surface. If this happens the 
position of the Fermi level at the surface shifts until the total charge 
on the surface, adsorbed ions, charge in the traps and charge in space- 
charge layer, adds up to zero. The consequences of this model have been 
carried through. 

The basis for the theory is illustrated in the energy level diagram of 
Fig. 13. At the semiconductor surface there is a space-charge layer of 
thickness Ca which gives a change in electrostatic potential of V B , 
corresponding to a potential energy of an electron of — eV B ■ Outside 
of the surface of the germanium proper, there is a surface film of thick- 
ness t D . A double layer giving a potential change V D , is formed from a 
charge of ions, 07 , on the outer surface of the film and charges in the 
surface traps of types a and b. Changes in c.p. with ambient result 
from changes in 07 ana! consequent changes in V B and V D . It is as- 
sumed that the remainder of the work function is independent of ambi- 
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Fig. 13 — Schematic of energy level diagram at germanium surface. 
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ent. When light shines on the surface, V B is changed to V B + &V B and 
there is an additional potential drop, bVi , in the body of the germanium 
resulting from the recombination current which flows to the interior. 
The change in contact potential with light is equal to SV B + 5F, . 

The film thickness ( D is shown on an exaggerated scale. We expect 
l D ~ 10~ cm and l B ~ 10~ cm, so that t B ^> l D . 

TABLE OF SYMBOLS 

A. Energies: 

E a = E a (true) + kT In (w u „ oc /co oc ), is the effective energy of the 
a-traps for V B = 0. Here co unoc and oj oc are the statistical 
weights of the unoccupied and occupied states, respectively. 

Eb = effective energy of the b-traps for V B = 0. 

E c = energy of lowest state of conduction band in interior of semi- 
conductor just beyond the space-charge layer. 

E v = energy of highest state of valence band at the same position. 

Ef = Fermi energy. 

Ei = E F when material is intrinsic. 

V D = potential drop across surface film. 

V B = potential drop across space-charge layer. 

Vno = value of V B for which n a = Pb , see below. 

Vq = value of V B o for an intrinsic sample. 

B. Concentrations: 

n = N c exp [(E F — E c )/kT] = equilibrium concentration 

(no./cm ) of conduction electrons in interior of semicon- 
ductor just beyond the space-charge layer. 

■p = N v exp [(E v — E F )/kT] = corresponding hole concentra- 

tion. 

n, = intrinsic concentration. 

n, , p a = equilibrium concentrations of electrons and holes, re- 
spectively, at the surface. 

N a , Nb = concentration (no./cm 2 ) of a- and 6-traps, respectively. 

n a = equilibrium concentration (no./cm ) of occupied a-traps. 

Pb = Nb — nb = equilibrium concentration (no./cm ) of unoc- 

cupied 6-traps. 

n a0 , Pbo = values of n a and pb for an intrinsic sample with V B = 0. 

ni = n + 5/i, pi = p -\- dp, and »« , p sl , n ai , p bi = concentra- 

tions in presence of light. Electrical neutrality requires that 
bn = 8p. 
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The theory is based on the following postulates : 

I. Changes in c.p. with ambient result from changes in 07 and the 
consequent changes in V B and V D . 

c.p. = V B +V D + const. (2) 

The charge 07 is largely compensated by charges in the surface traps. 
The barrier height, V B , is determined by the requirement of electrical 
neutrality: 

07 = e(n a — p b ) + const. (3) 

Since V B and V D are of the same order, the net charge per unit area 
in the space-charge layer will be smaller than 07 in the approximate 
ratio ( D /( B and may be neglected. 

II. Traps of type a have energies above E F and of type b below E F 
for all values of V B attained in the different ambients used. More ex- 
actly 

E a - eV B - E F > kT, (4) 

E F - E b + eV B > kT, (5) 

for all V B . One may then use the Boltzmann approximations for n a 

and pb '• 

Na 

Ua " 1 + exp [(E a - eV B - E F )/kT] (G) 

~ N a exp [(E F - E a + eV B )/kT\, 

Nb 

Vb ~ 1 + exp [{E F -E b + eV B )/kT] ( 7 ) 

~ N b exp [(E b - eV B - E F )/kT]. 

It is not necessary for our arguments to assume that all traps of each 
type have the same energy. The only requirement is that the distributions 
of trap energies are such that the a-traps are always above and the 6-traps 
always below the Fermi level for any ambient. 

III. Creation of electron-hole pairs by absorption of light occurs near 
the surface in a distance that is small compared with the diffusion 
length. Optical constants of germanium indicate that practically all of 
the light with energy sufficient to create electron-hole pairs is absorbed 
within a distance of 10 -4 cm of the surface. The diffusion length is of 
the order of 0.2 cm. 

IV. In the presence of light, the concentration of electrons in a-traps 
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is in equilibrium with the concentration of electrons in the conduction 
band and holes in b-traps are in equilibrium with holes in the valence 
band. The barrier height is adjusted so that the total charge in both 
types of traps is unchanged by illumination. We shall show in the ap- 
pendix that the resistance to flow of electrons from the conduction band 
across the space-charge layer and into a-traps is small compared with 
the resistance to flow of electrons from the valence band to the traps. 
Similar considerations apply to flow of holes to b-traps from the valence 
band as compared with flow from the conduction band. 

V. Recombination is limited by holes going into traps of type a 
and electrons going into traps of type b. The two types of traps act 
in parallel for recombination. The contributions to the surface re- 
combination velocity are proportional to p s n a and n,p b , respectively. 
These products are independent of V B and thus of ambient if postulates 
II and IV are satisfied. If other types of traps were important in re- 
combination, one would expect v, to depend on ambient, contrary to 
what is observed. 

Postulates I and II are used to relate V B and c.p. with changes in 
ambient. Postulates III, IV, and V are used to relate (Ac.p.) t with 
V B and the trap densities, and also to obtain an expression for the 
surface recombination velocity. 

As Vb is made more positive, corresponding to a decrease in barrier 
height for electrons, n„ increases and pb decreases. It will be convenient 
for the theoretical discussion to introduce the particular barrier po- 
tential V B o , for which n a = pb . With use of the Boltzmann approxi- 
mations in Equations (6) and (7), this gives 

exp [2/3 V B0 ] = (Nb/N a ) exp [(E a + E h - 2E,)/kT)] 

(8) 
= (Pbo/n a0 )(p/n), 

where /3 = e/kT. The last form follows from the definitions of pho and 
iiao and by noting that exp [2(E { — E F )/kT] = p/n. We then have 

njp b = exp [20(V B - V B0 )l (9) 

As so defined, V B o depends on the Fermi level and thus on the con- 
ductivity of the specimen. We shall let V be the value of V B o for an 
intrinsic specimen. From (8) 

n a0 /pbo = exp [-2/3Fo]. (10) 

If n a0 = pto , then V will be zero. Postulate II sets limits on V , but V 
is otherwise undetermined in our experiments. 
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We shall first relate V D and V B . Since the electric field in the film is 
A.iroi/K D , 

V D = 4*tna r /K D , (11) 

where K D is the dielectric constant of the film. Substituting for 07 
from Equation (3) and making use of (8) and (9), we find 

V D = 2H sinh 0(V B - Vbo) + const, (12) 

where 

H = *™k (N a N b exp [(E b - E a )/kT\) h = ^ire^/K D )(n a0 p b0 )\ (13) 

If V D is expressed in volts and l D in cm. 

ff = 1.8X 10" 6 (VJfi>)(»«opw)*. (13a) 

The contact potential Equation (2) may be expressed in the form 

c.p. = V B - Vbo + 2H sinh p(V B - V BO ) + const. (14) 

The change in contact potential with illumination results from a 
potential drop SVi in the interior and a drop 8V B across the space- 
charge layer. The former comes from the recombination current of 
holes and electrons diffusing from the surface to the interior. Since 
electrical neutrality requires that 8n = 8p, the concentration gradients 
are equal. However, the mobility of electrons is greater than that of 
holes, so that the diffusion current of electrons is larger than that of 
holes by the mobility ratio "b." Since there can be no net current flow 
to the interior, an electric field is established which is in such a direction 
as to enhance the flow of holes and retard the flow of electrons. The net 
potential drop associated with this electric field is 

SVi = 0>z« ln f, + (±±^\ (15) 

0(6 + 1) { bn + p j 

where 8p = 8n is the change in concentration in the interior just beyond 
the space-charge layer. 8 

This potential is positive for both n- and p-type material, and is 
independent of ambient, since 8p is independent of ambient. The ob- 
served (Ac.p.)z, is generally much larger than given by (15) and is of 
opposite sign for n- and p-type material. To account for the observa- 
tions, it is necessary to assume that the major part of the effect is 
associated with a space-charge layer at the free surface. We believe that 
the present experiments give the most convincing evidence obtained so 
far for the existence of such a space-charge layer. 1 
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The value of 8V B , the change in potential across the space-charge 
layer due to light, is determined by the requirement that there be no 
net change in charge in the surface traps, or that 

8n a = 8pb . (16) 

The changes 8n a and 8p b come both from changes in 8p and 8n and from 
8V B • According to postulate IV, n a \ = n a + 8n a is in equilibrium with 
the conduction band and p b \ — p b + &Pb with the valence band. 
We have then 

n a i/n a = n,i /n. = (ni/w) exp \08V B ], (17) 

Pbi/Pb = Psi/p. = (pi/p) exp [—ftiVa]. (18) 

from which it follows that 

&n a /n a = (ni/n) exp \06V B ] - 1, (19) 

8p b /p b = (pi/p) exp [-/9«7 B ] - 1, (20) 

* = exp MV B - V B0 )} = ^^^f 1 " 1 • (2D 
P6 (ni/n) exp [08 Kb] - 1 

Equation (21) is a quadratic equation in exp [/35F B ] which may be solved 
to give an explicit expression for 8V B • The role of electrons and holes 
may be interchanged by changing the sign of 8V B and of V B — V B0 . 
This accounts for the difference in behavior of n- and p-type samples. 
The total change with light is the sum of SVi and 8V B : 

(Ac.p.)z, = 8V L = 8V B + SVi . (22) 

In the analysis of the data, SVi is calculated theoretically from (15) 
with 8p = 8n determined from (1) and 8V B is obtained from the ob- 
served (Ac.p.)/, and SVi using (22). Equation (21) is then used to find 
V B — V B0 . A plot of c.p. — (V B — V B0 ) versus sinh (3(V B — V B0 ) should 
be a straight line with slope 2H. An analysis of the observed data in 
this manner is given in the following section. 

We turn finally to a discussion of the surface recombination velocity, 
v, . According to postulate V, recombination is limited by flow of holes 
to a-traps and of conduction electrons to fe-traps. The flow of holes 
from the valence band to a-traps (really electrons from a-traps drop 
into the vacant levels in the valence band corresponding to the holes) 
is proportional to the product of the hole concentration at the surface 
p,i , and the concentration of electrons in a-traps, n a \ . The reverse 
flow is that of thermal generation of holes: electrons from the valence 
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band go into unoccupied a-traps. Since, according to postulate II, the 
number of unoccupied traps is nearly equal to the total number of 
traps, and is thus independent of 8V B and 8p, the reverse flow will be 
practically equal to the thermal equilibrium value. If S pa is the re- 
combination cross-section for holes, the net flow of holes to a-traps is 

U pa = S pa v p (p a in a i - p.n a ) (holes/cm 2 sec). (23) 

Here, v p is the velocity factor which when multiplied by the concen- 
tration gives the number of holes crossing a unit area from one direction: 

v p = (kT/2irm p )K (24) 

With use of the relations, n ai /n a - n A /n, , p b \/pb = p,\/p* , from 
equations (17) and (18), which follow from the fact that the a-traps are 
in thermal equilibrium with the conduction band and the 6-traps with 
the valence band, (23) becomes 

U pa = S pa v P (n a /n,)(p a in,i — p,n,). (25) 

This expression may be simplified further. The ratio 



lln 



N a exp [(E F - E a - eV 



*vm , £■ exp ( h - F A (2fl) 
■ B )/m N c exp \ kT ) y } 



n, N c exp [(E F - E e - eV t 

is independent of E F and V B . The ratio may be evaluated for an intrinsic 
specimen with V B = 0, in which case n a = n a0 and n s = w< . Thus 

n a /n B = Uao/rii . (27) 
We also have from postulate IV, 

Psirisi = pini . (28) 
The equilibrium products are 

p,n„ = pn = Ui . (29) 
With use of (27), (28) and (29), (25) becomes 

U pa = S P aV p (nao/ni)(pini - pn). (30) 
Similarly, it is found that net flow of electrons to b-traps is equal to: 

U„ b = SnbVnipbo/ndipinx - pn). (31) 

The total rate of recombination is given by the sum of U nb and U pa 
and is given by an expression of the form : 

U = U nb + U pa = Cfanx - pn) = C(n + p)8p. (32) 
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It should be noted that the coefficient C is independent of the Fermi 
level and thus of the conductivity of the specimen, whereas v s is not. 
The relation between them is: 

For n-type C = v./n, (33) 

For p-type C = v,/p. (34) 

Values of C may be determined empirically from observed values of 
v 8 . Referring to Table I, for sample A, v a = 60 cm/sec and n = 1.4 X 
10 14 /cm 3 , so that C = 4.3 X 10~ 13 cm 4 /sec. For sample D, v s = 180 
cm/sec and p = 3.0 X 10 14 /cm 3 , so that C = 6 X 10 -13 cm 4 /sec. The 
values of C are approximately the same, indicating that the traps are 
not much different for the two specimens. 

The theoretical value of C involves Woo and pbo . It is the product 
n a pb = riaoPbo , which is related to the parameter H and which can be 
estimated from empirical data. To obtain the concentrations them- 
selves, a value must be assumed for V . We have 

n a0 = (noopw)* exp [— j87 ], (35) 

Pbo = (ttaopw) 1 exp [07 o ]. (36) 

As we have mentioned previously, there is no way to determine Vo 
from our experiments, although postulate II sets limits on its value. 

Let us for simplicity assume that S„bV n = S pa v p = S t v. Then, using 
(35) and (36) in (30) and (32), we have 

C = 2S«w(tt«oP»/tt*) i cosh PV . (37) 

This equation may be used to estimate the trapping cross-section S, 
for an assumed V . 

COMPARISON BETWEEN THEORY AND EXPERIMENT 

In Fig. 14 we have plotted c.p. — (c.p.) — (V B — V B o) versus 
2 sinh 0(Vb — V„o) for p-type sample D, (see Table I). Each symbol 
represents experimental points for one cycle. Results are shown for five 
different cycles not all in the same run. These results are typical of all 
the data obtained for this sample. It is seen that these data can be 
fitted over most of the range by a straight line as drawn, giving a value 
of H = 0.02 e.v. 

In Fig. 15 we have used the same experimental results to plot 
SV L versus c.p. — (c.p.) . Since the experimental values of 8V L cover a 
range of a factor of 20 or more, a logarithmic scale was used to plot 
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Fig. 14 — Plot of 2 sinh P(V B - Vi?o) versus c.p. - (V B 
as suggested by theory. 

8V L ■ Both positive and negative branches of the curve are plotted on 
the same figure. The symbols used for the experimental points are 
consistent with Fig. 14. When c.p. - (c.p.)o is greater than zero, 8V L is 
negative. As c.p. increases it approaches a maximum negative value, 
this is the negative branch. When c.p. - (c.p.)o is less than zero, 8V L is 
positive and as c.p. decreases 8V L approaches a positive maximum 
value that is less in magnitude than the negative maximum. The solid 
curve represents the prediction of theory for H = 0.02 e.v. The agree- 
ment between theory and experiment is good. It should be emphasized 
that this fit is obtained with only one adjustable parameter. 

The data for the other samples were analyzed in the same way. 
The results are shown by plotting 8V L versus c.p. - (c.p.)o as in Fig. 15. 
Fig. 16 is for n-type sample A and Fig. 17 for w-type sample B. The 
values obtained for H were 0.015 and 0.022 e.v. respectively. The fits 
obtained are about equally good in all cases with some deviation between 
theory and experiment near the extremes of contact potential. The 
values of H obtained are all of the same order as they should be if the 
surface trap structure is approximately the same from sample to sample. 
When V B - V B o is large and positive and thus c.p., Eq. (14), is large, 
8V L approaches the negative maximum 
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In — 



)-«]■ 



equation 21, and as V B — V B o becomes small and negative, c.p. decreases 
and bV L goes through zero and approaches the positive maximum 



+ 



[G-?) + *4 



For p-type germanium njn is greater than pi/p so that the negative 
maximum for 8V L is greater than the positive maximum. Just the 
opposite is true for ?i-type germanium. 

In this comparison of theory and experiment a tacit assumption has 
been made that the germanium surface is uniform in its properties. It 
might well be that this is not the case. The surface might be "patchy." 
To estimate what effect patches or non-uniformities might have the 
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Fig. 15 — Change of contact potential with illumination, iVj. , versus c.p.; 
experiment and theory sample D, p-type. 
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Fig. 16 — Same as Fig. 15 for sample A, n-type. 

following was done. It was assumed that the surface was made up of 
two parts each having the theoretical dependence of 8V L on c.p. but 
with the contact potential where the light effect goes through zero 
differing by 0.05 volts. The values of 5V L at a given contact potential 
were averaged and plotted against c.p. The difference between this 
curve and the original one was almost entirely a simple shift in (c.p.)o • 
It is therefore unlikely that non-uniformities in c.p., over the surface, of 
the order of 0.1 volts or less would change the relation between 8V L 
and c.p. sufficiently to be detectable. 

One can use equations (8), (10) and (21) to calculate values of 
(y B _ v ) for each experimental value of 8V L and then plot these 
values against the corresponding values of contact potential. These 
results are shown in Fig. 18 for w-type sample A and p-type sample D. 
(V B — V ) changes with c.p. as one would expect and in an approxi- 
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mately linear fashion throughout most of the experimental range. The 
change in V B is approximately one-fifth the change in contact potential. 
(V B — Vq) is positive for the p-type sample and negative for the n-type 
sample throughout most of the range. If the trap distribution on the 
surface were symmetrical about the intrinsic position of the Fermi 
level, then pw would be equal to n a0 and V would be zero. In this case 
the space charge layers on n- and p-type germanium would be about 
equally developed. All the experimental evidence on germanium indi- 
cates that this is not the case but rather that — V Bn is much larger than 
V Bp ■ This then means that Vo is negative and that the trap distribution 
is unsymmetrical either in number or energy or both in such a way that 

>>a0 > Pbo . 

In Table II values are given for the differences (V B — V B o), etc., 
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Fig. 17 — Same as Fig. 15 for sample B, n-type. 
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-0.06 



0.48 0.52 0.56 0.60 0.64 0.68 
CONTACT POTENTIAL IN VOLTS 



0.72 



0.36 0.40 0.44 

Fig. 18 — Potential across space charge layer, V B , versus contact potential. 



0.80 



also 8V's, for certain special cases for samples A and D as calculated 
from theory. From equations (8) and (10) 



(J bo - V ) = - (E> - Er) = ±- ln(p/n). 
e £p 



(38) 



It is a constant for any given sample as is 8Vi (see equation 15). From 
equations (21) and (38) it follows that when 8V B = then (V B - V ) = 
for all samples. To say this another way, if the traps are symmetrically 
distributed about Ei and V B is zero, then illuminating the germanium 
would produce no potential difference across the surface. One would not 
suspect offhand that (V B — V ) for 8V L = was also approximately 
independent of the material in the sample but it is. For the case where 
8p and the minority carrier are both small compared to the majority 

Table II 





A (n-typc) 


D (f-type) 


(Vbo - Vo) 


-0.044 


+0.064 


SVi 


+0.0019 


+0.0015 


{V B - V a )&v B = 








(V B - V B0 )&v B = 


+0.044 


-0.064 


(V B - V )&v L = 


+0.010 


+0.010 


(V B - V B0 )iv L = 


+0.054 


-0.054 


SV L , (V B - V BO ) = 


+0.024 


-0.026 
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carrier one can easily show that when 8 Vz, = 0, 
(V B - V ) = (1/2/3) ln6. 

This follows from equations (15), (21), (22) and (38). The rest of the 
results hardly need comment. 

It is interesting to compare differences in contact potential, etc., for 
units A and D. In equation (2) the Fermi energy has been included in 
the constant. This equation predicts that if the trap distributions on 
both surfaces are the same, then when the contact potentials of both 
surfaces are equal the values of V D should also be equal and the dif- 
ference between the TVs should be equal and opposite to the differences 
between the Fermi energies in electron volts. This equation can be 
written 

c.p. = (Vb - V ) - (V B0 - Vo) + V D 4- V + const . (40) 

In comparing units we shall use A's to denote differences and these 
differences are always taken A-D. 

For the case where Ac.p. = 0, A(V B - V ) can be read from Fig. 18 
and A(V B o - V ) from Table II. We have then 

A(V D 4- V + const) = -0.05 . 

This indicates that our simple picture is not quite right. If AV D were 
zero for this case it should follow from equation (12) that both 
A2H sinh /3(V fl — V B o) and the difference in the const in this equation 
should be zero. It turns out, however, that A2H sinh /3(F B - V B o) is 
not zero but +0.07. Equation (12) can be substituted into equation 
(40) giving 

c.p. = (V B - V ) - (Vbo - Vo) + 2H sinh (*(V B - V BO ) 

+ V + const 2 . 

Where the constant now includes the constant part of V D and is labeled 
const 2 to distinguish it from the constant in equation (40). We have 
then for Ac.p. = 

A(7 + const 2 ) = -0.12. 

This difference A(V 4- const 2 ) can be calculated three other ways, 
using the experiment results and theoretical values where necessary. 
These ways are 
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1) at the same time in the cycle 

A(V + const 2 ) = -0.10 

2) when 8V L = 

A(V + const 2 ) = -0.13 

3) when 8V B = 

A(V + const 2 ) = -0.12 . 

All four results are consistent within the probable experimental accuracy 
and give an average result of A(V + const*) = —0.12 instead of zero 
as one might expect from the simple picture. This indicates that the 
trap distributions are different on the two surfaces. The constant part 
of V D is proportional to — (N a — N b ) and from equation (10) one would 
expect V to decrease as N a increases with respect to N b , thus AV 
and A const 2 should be of the same sign and additive so that A const 2 
is less than 0.12 volts; assuming that ( D /K D is the order of 2 X 10 -7 cm, 
this indicates that A(N a - N b ) is the order of 3 X 10 11 per cm 2 which is 
small compared with the probable trap density, N a , N d , of the order 
of 10 M per cm 2 as we shall see in the next paragraph. 

Assuming that C D /K D in equation (13a) is the order of 2 X 10" cm 
one can calculate (rupw))* obtaining 4.1 X 10 10 and 5.5 X 10 for 
samples A and D respectively. Using these values one can solve for S t 
in equation (37) and obtain for the case of V = the value 5.0 X 10~ 17 
cm 2 for the average capture cross section of the surface traps. The 
values of S t , n al) and p M depend on what one takes for V . The relations 
are 

5 X 10~ 17 
1 cosh Wo ' 

n aQ = 5 X 10 ,0 exp[-/3Fo], 

Pb0 = 5 X 10 10 exp \fiVo] • 

This dependence is shown in the graph in Fig. 19. As already mentioned 
there are reasons for thinking that Vo is less than zero. If one takes 
— 0.06 volts as a reasonable value then one gets S t =10 cm , riot, 
= 5 X 10 u /cm 2 and p a0 = 5 X 10 9 /cm 2 respectively. One can push these 
calculations still further to estimate N a , N b , E a and E b . One knows 
that E a and -E b must be greater than 1/kT. Also JV« and #6 should be 
less than the number of germanium atoms per cm 2 of surface which is 
1.4 X 10 15 /cm 2 . Values of N a and N b of the order of 1 X 10 14 /cm 2 for the 
number of traps per cm 2 with energies E a and —E b of the order of 0.2 
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e.v. measured from the midband energy are reasonable and not incon- 
sistent with the original assumptions. 

As mentioned in the experimental section, some data on (Ac. p.)*. for 
samples C and E have been obtained. While no complete analysis of 
these results has been made, one can see that they are of the right 
order of magnitude. As the specific resistance of the sample decreases, 
the body life time t decreases. This empirical result is to be expected. 
Consequently 8p for the same light intensity decreases and one would 
expect (Ac.p.)i, to decrease as it does. For sample E of course one could 
not neglect the charge in the space layer so that the theory would be 
more involved. 

The comparison between the contact potentials of samples A, C, D 
and E shown in Fig. 6 can be understood in part at least. Consider first 
the over-all result that at the same time in the cycle the c.p. for a sand- 
blasted surface is less than for an etched surface, i.e., the work function 
for the sand-blasted surface is greater. It is known that the surface re- 
combination increases enormously when the surface is sand-blasted. This 
means that either the surface trap density has increased or that the 
distribution has changed or both in such a way as to increase surface 




-0.08 -0.06 -0.04 



-0.02 0.02 

Vn IN VOLTS 



0.04 0.06 



Fig. 19 — Dependence of currier densities in surface traps, n a0 and Pbu , for 
Vb = 0, and cross section of trapping S t on V . 
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recombination. The results of Fig. 6 indicate that this change in trap 
distribution also has been in a direction to increase (— V B ). This is what 
one would expect if N b has increased with respect to N a , i.e., sand-blasting 
increases the effective number of acceptor traps. About all that can be said 
about the rest of the results in Fig. 6 is that if one knew the surface trap 
structure in number and energy distribution for the etched surface then 
one could deduce from the results in Fig. 6 the distribution of traps for the 
sand-blasted surface over part of the energy range at least. 

The theory developed here for the germanium surface may not apply 
at all for a silicon surface. In a previous discussion of the silicon surface 
it was assumed that the resistance to surface trapping occurred mainly 
in the flow across the space charge layer rather than from the surface 
to the traps. This may well be the case in silicon. An investigation of the 
silicon surface using these same techniques should clarify this point. 

It should be emphasized that the ambient used for this study and the 
variations in it are special in that they do not necessarily correspond to 
the atmosphere of ordinary room air. It is very probable that constituents 
in room air other than oxygen ions and water vapor are important, such 
as salt ions for instance. It is quite probable that v, for such a surface 
exposed to room air would increase considerably with time instead of 
staying relatively constant as it does under the bell jar. 

CONCLUSIONS 

A method has been developed for studying the surface properties of 
Ge in a gaseous ambient at atmospheric pressure. It has been found 
that the Ge surface interacts with this ambient. Two atmospheric 
constituents that are important in this interaction, oxygen and vapors 
with OH radicals, have been isolated. With the controlled use of these, 
the surface dipole of Ge can be cycled between two extremes. Thus the 
dependence of other properties of the Ge surface on surface dipole, such 
as change of contact potential with illumination and surface recombina- 
tion, can be determined. 

It is evident from the results that the method is very powerful. In 
order to complete the present study, it was necessary to stop trying new 
experiments since almost all directions of variation open up new and 
interesting phenomena. These results on Ge are really just a beginning, 
and the preliminary data on silicon indicate that the same method would 
also be fruitful on other semiconductors. The technique is of course not 
limited to atmospheric pressure. 

A tentative theory of the Ge surface has been developed that is 
sufficient to explain the experimental results on a semi-quantitative 
basis. Theory and the experiment together predict approximately the 



SURFACE PROPERTIES OF GERMANIUM 37 

number, type and distribution in energy of the surface traps. One has 
therefore a tentative model of the Ge surface that should be very useful 
in any further investigation of its properties. 
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Appendix 

We have assumed (Postulate IV) that traps of type a are in equi- 
librium with electrons in the conduction band and that traps of type b 
are in equilibrium with holes in the valence band. We wish to show that 
this is not really a separate assumption but follows as a consequence of 
Postulate II if the density of traps is not too high. For simplicity, we 
shall restrict the discussion in the appendix to the limiting case of small 
departures from equilibrium so that the equations are linear. The 
problem may then be discussed most conveniently by means of quasi- 
Fermi levels* or imrefs, 4> n and <f> p , for the conduction electrons and 
holes, respectively. 

Departures 5<£„ and 84> p , of the imrefs from the Fermi level are a 
measure of the departures of the concentrations, 8n and 8p , from their 
equilibrium values: 

fit = n + &n = n exp (— 08<f>„), (A.l) 

Vi = V + 8 P = V exp (Pfyp). (A.2) 

The imrefs in the interior are then 

* - T '" ( x + t) • (A ' 3) 

8 *> - \ h + 7) • (A ' 4) 

Correspondingly, the imrefs at the surface are defined by: 

n A = n, + 8n, = n exp [0(V B + 8V B - M> ns )], (A.5) 

p* = p, + 8p, = v exp [-/3(7 fl + 8V B - ftfe,)]. (A.6) 



* Reference 2, pages 302-308. 
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Changes in imrefs of the traps are defined by: 

nn = n t + 8n t - j + ^ [( ^ _ & _ e( ^ + ^ + e ^ )/A . r] . ^-" 

pa = pi + ty« = i + exp [( Ef - I t + e(V B + 5F fl ) - e8<f> t )/kT] ' (A>8) 

In these last two equations, t may refer to either type of trap (a or b) 
and 8p t = —8n t . 

For small signals, the recombination current via a given set of surface 
traps may be considered as a flow produced by differences in the imrefs, 
&/>„ and 84> p , through four effective resistances in series, as shown in 
Fig. 20. Here R nB is an effective resistance for flow of electrons across 
the barrier layer, R„ t for flow of electrons from the conduction band at 
the surface to the traps, R pt for flow of holes from the filled band to the 
traps, and R pB for flow of holes across the barrier layer. Under steady 
state conditions, the net flow of conduction electrons to the surface is 
balanced by an equal flow of holes. The recombination current may be 
thought of as a flow of electrons from the conduction band via the traps 
to the valence band. 

The recombination current per unit area is 

/ = -eU = (8<t> n - 8<t> p )/R t , (A.9) 

where R t is the sum of the four resistances in series. Here U is the 
particle current and I is the corresponding electric current. If 8<t> n and 
8<t> p are expressed in volts, R t is in ohms/cm 2 . We may define a re- 
combination constant Ct by an equation corresponding to (34) : 

U = Ct(pmi - nd. (A. 10) 



6f t 




Fig. 20 — Circuit analogy of surface recombination. 
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The relation between C t and R, is obtained by use of (A.l) and (A.2) 
in (A. 10). Since (A.9) is valid only to terms of the first order in ty„ 
and 8<f) p , we make the corresponding approximation in (A. 10) and find 

1/7?, = ePn*C t . (A.11) 

If there is more than one type of trap, the net recombination resistance, 
R, is that of the various traps in parallel. Values obtained for R for 
specimens A and D from the empirical values of C are about 500 ohms. 

We shall show that for a-traps, R pa is much larger than the other 
resistances in series, and that for o-traps, R nb is the dominant resistance. 
This implies that 60 a = 50„ and 8<t> b = 50 p , or in other words that a-traps 
are in equilibrium with the conduction band and b-traps with the 
valence band. 

First consider the flow across the space-charge layer. The resistances 
R nt and R bl depend on the sign of V B . When V B is negative, the net 
electron current across the space-charge layer is: 

/ = ev n (n al - n) exp \fi(V B + &V B )] 

^ ePv„n a (8<t> n — 50 ns ), 

where v n is defined by an equation similar to (24). The second form is 
the linear approximation. Thus we have 

1/RnB = e0v n n. . (A.13) 

If V B is positive, n, is replaced by n. 

Correspondingly, for V B negative, the hole current is: 

/ = cv p (pi - p A ) exp \fi(V B + 57b)] 

^ epv p p(84> ps - 8<t> p ), 

and 

1/R p b = e$v p p. (A.14) 

For V B positive, p is replaced by p s . 

The maximum value of R nH is obtained for the ambient which makes 
V B most negative and R pB is a maximum when V B is most positive. 

The current from the conduction band to the traps is calculated as in 
(23): 

/ = e(v n S n tn tl p,i - g t n n ) (A. 15) 

£ eBv n S n fit.p t (84,n. - 60i). (A.16) 
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Here g t is the rate of thermal generation of electrons from occupied 
traps to the conduction band. Since, for equilibrium conditions, 1 = 0, 
we have: 

gt = v„S„tn a pt/n t . (A.17) 

It may be verified easily that g t is independent of E F and V B . The 
ratio, 

n ct = n s p t /n t = N c exp [(E t - E c )/kT], (A.18) 

is the equilibrium concentration of electrons in the conduction band 
when the Fermi level is at the level of the traps. The expression (A. 16) 
is again the linear approximation. Thus 

l/R K t = ePv n S nt n s p t . (A.19) 

Similarly, we have 

1/R p t = e(3v p S pt p,n t . (A.20) 

We shall show" that the ratio 

Rm/Rpt = (vnS p t/v p Snt)(p,nt/n,p t ), (A.21) 

may be expected to be small compared with unity for a-traps and large 
compared with unity for 6-traps. First consider a-traps. If anything, 
S pa < S na , because holes must give up a larger energy than conduction 
electrons in going to a-traps. The second factor is small compared with 
unity if 

p s « n ct , (A.22) 

with net defined by (A.18). This will be the case if the a-traps are closer 
to the conduction band than the top of the valence band at the surface 
is to the Fermi level. According to Postulate II, this should always be 
the case. 

Similarly, for 6-traps 

if 
where 



Rpt/Rnt « 1 (A.23) 

n, « p vt , (A.24) 



p v t = p.n t /p t = N v exp [-(E t - E v )/kT]. (A.25) 

is the concentration of holes in the valence band when the Fermi level 
is at the level of the traps. 
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The values of R„ B and R pB relative to R nt and R pt are: 

RnB/Rnt = SntPt (A.26) 

R pB /R P t = S pt n t . (A.27) 

In our experiments these are always small compared with unity, since 
the trapping cross-sections are of the order of 10 -17 cm 2 and p t and n t 
are always less than Nt , which is of the order of 10 14 /cm 2 . 

Barrier resistances may be important for surfaces with a large number 
of surface traps. Analysis of earlier data on the change of the contact 
potential of silicon with light 1 was made on the assumption that the 
probability that an electron or hole reaching the surface be trapped is 
relatively large, and that consequently the barrier resistances are large 
compared with the trapping resistances. The present experiments on 
germanium throw doubt on this interpretation, but further experiments 
are required to clarify the situation. 
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